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Abstract 

A uniform minimal model of rhodopsin photoisomerization induced by either coherent 
laser light or low level incoherent light (e.g. moonlight) is provided. Realistic timescales for 
both processes, which differ by ten orders of magnitude, are obtained. Further, a kinetic 
scheme involving rates for both coherent and incoherent light excitation is introduced, plac- 
ing all timescales into a uniform framework. 
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I. INTRODUCTION 

Developments in fast pulsed lasers have allowed for the detailed study of photobiological 
processes such as laser induced cis /trans isomerization of rhodopsin, a process of interest 
due to the large quantum yield (~ 65 %), high speed (~ 200 fs) of reaction, and importance 
in the function of living organisms l|, |2|, |3] • However, photoinduced processes such as this 
occur naturally in the presence of weak incoherent light, rather than in the strong coherent 
light that emanates from laser sources. For example, photoabsorption in rhodopsin initiates 
vertebrate visual transduction in dim light, such as moonlight [4j. Since the processes 
induced by these two types of sources are qualitatively different, e.g. pulsed coherent light 
induces time dependent molecular dynamics, whereas purely incoherent light does not 5j, |6j , 
it is important to establish the relationship between them. 

In this paper, we provide a uniform minimal model for photoisomerization induced by 
either of these light sources and demonstrate: (a) a computed dynamics timescale for fem- 
tosecond laser pulse excitation in agreement with experiment, (b) realistic dynamics for 
time scales on the order of milliseconds for moonlight induced processes, and (c) a kinetic 
scheme involving rates of both incoherent and coherent excitation that places all timescales 
within a unified framework. Specifically, in the natural visual process, the femtosecond co- 
herent timescales provide the initial rise of the cis/trans isomerization and the millisecond 
incoherent timescale gives the rate of the process at longer times. 



4 



II. THEORY 



Our theoretical treatment of the photoisomerization is based on a one dimensional system 
with two electronic states (see Fig. UK) connecting the cis and trans configurations, coupled 
through a strength parameter 77 to a "bath" that models the effects of the remaining degrees 
of freedom and of the external environment. Isomerization occurs via rotation about an 
angle a. The interaction potential between the system and the coherent external field 
E(t) is treated by means of the dipole approximation. In the case of low level incoherent 
light, E(t) = and a second bath describing the incoherent light is included. That is, our 
Hamiltonian is 



where H$ is system Hamiltonian, /i is transition dipole moment of the system, E(t) is 
electric field of the laser pulse, H env is the environment Hamiltonian, i^ienv is the interaction 
Hamiltonian between the system and environment, if ra d describes blackbody radiation, and 
ifi ra d is interaction Hamiltonian between the system and the radiation field. Eigenstates \i) 
of the system H$ satisfy 



and the density matrix accounted with evolution of the (system + bath) is denoted px- The 
system density matrix is p = Ty-qPt, where Ttb denotes a trace over the bath. The time 
propagation of the density matrix elements of the system Pijit) = (i\ p(t) \j) is described by 



H T = H S — pE{t) + iflenv + H env + Hi rad + -f/rad, 



(1) 



H s \i) = Xi \i) 



(2) 



Redfield theory within a secular approximation [7|, |8|, |9|, [10j as 
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where Wji = T^,- + 1^ . ^ is transition probability per unit time from ith to jth eigen state 
of H s , and jij = J2 k i T tkki + T jkki) ~ T tju ~ T m is dephasing rate. Here, 
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where the brackets (. . . ) B represent a trace over degrees of freedom in B, where B is either the 
environment "env" or the incoherent radiation field "rad", and Hi^(t) = e lHBt ^ h Hi^e~ iHBt ^ h . 



The system Hamiltonian Hs is given in terms of two diabatic electronic states by 

fT + V g (a) V ge (a) \ 
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is the kinetic energy, V g (a) and V e (a) are the potential energy surfaces 



in ground and excited electronic state, and V^ e (a) = V cg (a) is the coupling potential between 



ground and excited states (see Fig. [T^l) . 



The environment is described as a set of harmonic oscillators of frequency uo' n and the 
system-environment coupling is H\ eny = QJ2 n ^ K n(bn + b n ), where tf n and b n are the 
creation and annihilation operators pertaining to the nth harmonic oscillator. The op- 
erator Q is a diagonal 2x2 matrix with cos a on the diagonal, and the coupling con- 
stants K n and spectrum of the bath are chosen in accord with an Ohmic spectral density 
J(uj) = 2ir J2 n K n!^( u} ~ w n) = V 1 ^ 6 -^ i where the strength of the system-environment cou- 
pling is determined by the dimensionless parameter rj, and u c = 300 cm -1 . After some 
algebra, we obtain first term of Eq. (jSJ) as, 
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where n(u) = {exp(faj / k^T) — 1}~ is the Bose distribution at temperature T = 300 K , 
Uji = (Xj — Aj) /h, and Aj is an eigenenergy of H s . 

As a typical situation of scotopic vision, we consider moonlight, which is well characterized 
as a blackbody source at 4100 K [HI]. The radiation field is also described as a set of harmonic 
oscillators of frequency u;" and the system-radiation field coupling is treated by means of 
dipole approximation as, 



trio" 

#irad = 2eV sin61 ( afcex P ( ik ' r ) ~ flfct ex P (~ ik ' r )} ' ( 9 ) 

where k is a wave number vector, eo is the permittivity of vacuum, r is a position inside of 



a cavity, V is vo 



ume of the cavity, and 9 is an angle between the transition dipole moment 



vector and k 12j. By assuming the large cavity limit the summation of k can be replaced 



with integrals, and second term of Eq. (jSJ) is written as, 
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A component of the imaginary part of Eq. (1101) describes the Lamb shift. The integration 



with respect to k does not converge, and this difficulty can be avoided by renormalization 



theory 



131 ]. However, since the effect of Lamb shift is generally less than 0.1cm 1 , the 



divergent term in Eq. ffTO]) is neglected in this paper. The coefficient C in Eq. fflUj) is 



introduced to adjust density of blackbody radiation to that of light incident on our retina. 
Specifically, by assuming that one is looking at a surface lit by moonlight, with a color 
temperature of 4100 K and a luminance LCd-m~ 2 , the ratio of the intensity of light falling 
on the retina over the light falling on the cornea as 0.5, the pupil area 3.8 x 10~ 5 m 2 , and 
the distance from the lens to the retina of 0.0167 m, we obtain C = L/4.0 x 10 10 . Here, 



a conversion from luminous flux in Cd-sr to radiant flux in W-m 1 was done by using the 



spectral luminous efficiency function for scotopic vision 14|. 



From Eqs. (|SJ) and ffTO]) . we obtain the transition probability in (jSJ) as, 
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where A+j and Bij are Einstein A and B coefficient in between the ith and jth eigenstate of 
H$, and W{oS) is the Planck's energy density. The dephasing ratio 7^ in (j3J) is evaluated by 
numerical integration of Eq. (JHJ) . 

III. RESULTS AND DISCUSSION 

Figure [lb shows the time propagation of molecular populations under a typical laser 
pulse of time duration 5fs, amplitude 4 x 10 9 V/m, and a carrier frequency of 2 x 10 4 cm _1 
that is resonant with the excitation to the electronic excited state around the Franck- 
Condon region. The transition dipole moment, set at lODebye, corresponds to an oscillator 
strength / rs 1. At time t = 0, the cis population P C i S (t = 0) is almost unity and after 
t = 10 fs, probability is created in the excited state. Each panel in the Fig. [Tb shows the 
relaxation process with a different degree of system-environment coupling: rj = 12.5, 25 and 
50. Evident is the fact that the trans yield is lower, and the isomerization is faster, with 
increasing coupling rj to the bath. We note that the time scale of the reaction in Fig. [TJ is in 
accord with that observed experimentally using coherent light excitation of rhodopsin, i.e. 
on the order of 200 fs (l|,b|. 

By contrast, the time dependence of the molecular populations for the case of excitation 
by incoherent light is shown in Fig. [2J Here we examine the problem in a context relevant 
to realistic biological systems. As seen in Fig. [2], for all 77 the rate of increase of Ptrans is 
linear in time after a time that we denote as t c {rf). Subsequent to that time the slope of 
Ptrans vs. t is s — 9.4 x 10~ 8 s~ 1 , corresponding to a cis/trans isomerization timescale of 



almost one year. Note that the slope s is independent of the speed of photoisomerization 
observed under pulsed laser conditions, as evidenced by the fact that it is independent of rj. 
Rather, this rate of transformation is dictated by the photon flux, which is the rate limiting 
reagent in the process. By contrast, the time t c , which corresponds well to the time scale of 
photoisomerization under the laser pulse, relates directly to rj as t c r\ ~ 20 ps. For example, 
for the case of r\ = 12.5, t c = 1.5 ps, in accord with Fig. [lb. 

Figure [2b shows the time dependence of Ptrans as a function of the luminance L of the 
incoherent light source. The slope s is seen to be proportional to the luminance L as 
s/L « 3.1 x lO^Cd-W-s- 1 . 

Since the isomerization of only a few molecules are necessary to induce hyperpolarization 
in a rod cell [4], we compute Ps, the probability that at least three from among all of the 
cis molecules in a rod cell are converted to trans. The probability would then correspond 
to the rate of our initial visual process under moonlight conditions. The probability -P^t) 
that at time t at least three from among N molecules are trans is given by 1 — po — p\ — P2, 
where 

p n = C n N p n {l-p) N - n (12) 

is a probability that n from among N molecules are converted to trans. Here, p = Ptranlt) is 
the probability that a molecule is trans at time t, and is the binomial coefficient. For the 
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case of vision, we take the number of rhodopsin molecules in a rod cell to be iV = 4x 10 
and assume that the time dependence of Ptrans maintains a constant slope s until t = 25 msec. 
The resultant P3 values are shown in Fig. [31 where the time scale to obtain at least three 
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trans molecules is on the order of a few tens of milliseconds. This finding is consistent with 
experimental time scales of 10 msec for dim flash response of a rod cell 4|. We note, as 
in the previous results, that the speed of photoisomerization under pulsed laser conditions 
bears no relation to the far longer time scales associated with the evolution of probability 
P3, since the photon flux is rate-determining in the latter case. Note further that the times 
at which Ps(t) reaches the value of 0.5, a measure of the biological response, is virtually a 
linear function of the irradiance. 

Thus far, molecular time evolution in incoherent light was considered using the Redfield 
approach. We also find that the population transfer can be modeled analytically by solving 
the simple three state model with the four reaction rates shown in Fig. HI A comparison with 
the computed Fig. [2] gives excellent results. Here, states A, B, and C represent cis, excited, 
and trans conformations of the molecule, respectively. The values of k 2 and k^ correspond 
to rates of population transfer from P c to P C i S and Ptrans, which are mainly caused by the 
system-environment coupling. Values obtained from the coherent pulse studies of Fig. [1] 
give k2 = & 4 = 0.08r7 ps" 1 . The ki and k% represent rates of population transfer from 
P cis and Ptrans to P e , caused by both system-environment coupling and photoabsorption. 
The rates of system-environment coupling can be assigned using detailed balance, and the 
rates of photoabsorption are given by the Einstein transition probability from the electronic 
ground state to the electronic excited state. In the case of k\, the primary contribution is 
photoabsorption, giving k% = BW = L x 5.6 x 10~ 6 Cd _1 -m 2 -s _1 , where B is the Einstein 
B coefficient, and W is density of energy of the radiation field. The densities of the field 
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used in Fig. @] correspond to the luminescence values used in Fig. [2] 16]. On the other 
hand, in the case of k^, the dominant term is system-environment coupling, and we obtain 
&3 = &4 x 1.87 x 1CT 9 . With the resultant ki, k 3 « k 2 , k±, the rate equations give 
the reaction rate for isomerization under incoherent light as ki/2 = BW/2. Further, these 
equations establish the existence of a linear region for Ptmns vs. t with an r\ independent slope 
s = k\/2 after a time t c = 3/(k 2 + k^), relating the rate approach to both the computed 
coherent and incoherent results. 

We note that the reaction rate obtained by the three state model is ~10 % smaller than 
that given by the Redfield equation. The difference mainly comes from the simplifying 
assumption that k 2 = k±, and the evaluation of the rate of photoabsorption at the torsional 
angle a set to zero. Nonetheless, all of the trends seen in the Redfield computed results are 
also evident in the rate equation results. 

IV. SUMMARY 

We have presented a unified theoretical model of photoisomerization under both a co- 
herent light source such as a femtosecond laser pulse and an incoherent light source such as 
the moonlight. A minimal model of the isomerization process that gives the same timescale 
as the femtosecond laser experiment was obtained. It was shown that the time scale for 
photoisomerization under coherent light corresponds to the initial rise time t c of the photoi- 
somerization under incoherent light. Further, we introduced a simple three state model that 
incorporates all of the relevant rates obtained from both the femtosecond and millisecond 
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time domains. 

This approach provides a connection between the time domain of the femtosecond laser 
experiment and that of biologically relevant response time scales. A dynamical behavior 
is seen even for the case of incoherent light source, since sudden irradiation of the light at 
t — introduces partial coherence into the system. The very earliest dynamics correlate 
with the primary event of isomerization as identified in femtosecond laser experiments. The 
exact response observed reflects the combined effect of the characteristics of the radiation 
field and the underlying dynamics. 
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Figure Captions 

FIG. [1] a) Potential energy surfaces for the two state model for cis to trans photoisomer- 
ization. The solid curve and dotted curve show diabatic potentials and V e , respectively. 
The dashed curve shows a coupling potential between two diabatic electronic states, b) Time 
propagation of cis and trans populations under a short intense pulse for different values of 
r\. P cis is the population in the range — | < a < | on V g , F \ rans is that in the range 
— 7r < a < — on V e , and P c = 1 — P C i S — P trans . Note that in Panels b and c, the very short 
time dynamics, which includes the excitation from the cis, is not evident due to the short 
time over which it occurs. 

FIG. HJ a) Time dependence of Ptrans for three r\ values with incoherent light luminescence 
L = 0.03 Cd-m -2 . b) Time dependence of Ptrans for various values of L, with system- 
environment coupling 7] = 25. In all cases there is a deviation from strictly linear behavior 
at the early times that corresponds to timescales of isomerization dynamics. 

FIG.|3j Time dependence of the probability P3 that at least three from among 4xl0 9 cis 
molecules become trans for various values of L: 0.060 Cd/m 2 (solid); 0.030 Cd/m 2 (dotted); 
0.015 Cd/m 2 (dashed). 

FIG. SJ a) Three states model with reaction rates k±, k 2 , k 3 , and k^. Time propagation 
of Pc for each luminance L and strength parameter of system-environment coupling rj. 
Compare with results shown in Fig. [2j 
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